A Generic Dynamical Model of Gamma-ray Burst Remnants 
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The conventional generic model is deemed to explain the dynamics of 7-ray burst remnants 



adiabatic expansion, the model could not reproduce the Sedov solution in the non-relativistic 



phase, thus the model needs to be revised. In the present paper, a new differential equation 
is derived. The generic model based on this equation has been shown to be correct for both 



radiative and adiabatic fireballs, and in both ultra-relativistic and non-relativistic phase. 
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1 INTRODUCTION 

Since the BeppoSAX detection of GRB 970228, X-ray afterglows have been observed from about 
15 gamma-ray bursts (GRBs), of which 10 were detected optically and 5 even also in radio 
wavelengths (Costa et al. 1997; Kulkarni et al. 1998; Bloom et al. 1998; Piran 1998; and 
references therein). The cosmological origin of at least some GRBs is thus firmly established. 
The so called fireball model (Goodman 1986; Paczyhski 1986; Rees & Meszaros 1992, 1994; 
Meszaros & Rees 1992; Katz 1994; Sari, Narayan & Piran 1996) is strongly favoured, which 
is found successful at explaining the major features of the low energy light curves (Meszaros 
& Rees 1997; Vietri 1997; Tavani 1997; Waxman 1997; Wijers, Rees & Meszaros 1997; Sari 
1997; Huang et al. 1998; Dai k, Lu 1998a; Dai, Huang & Lu 1998). A variant of this model, 
where central engines (e.g., strongly magnetized millisecond pulsars) supply energy to postburst 
fireballs through magnetic dipole radiation, has been proposed to account for the special features 
of the optical afterglows from GRB 970228 and GRB 970508 (Dai & Lu 1998b, c). 

Since the expansion of a fireball may be either adiabatic or highly radiative, extensive 
attempts have been made to find a common model applicable for both cases (Blandford & 
McKee 1976; Chiang & Dermer 1998; Piran 1998). As a result, a conventional model was 
suggested by various authors (see for example, Chiang & Dermer 1998; Piran 1998). A dynamical 
model should be correct not only in the initial ultra-relativistic phase, which is well described 
by those simple scaling laws (Meszaros & Rees 1997; Vietri 1997; Waxman 1997), but also in 
the consequent non-relativistic phase, which is correctly discussed by using the Sedov solution 
(Sedov 1969; Wijers, Rees & Meszaros 1997). Although the conventional model is correct for 
the ultra-relativistic phase, we find it could not match the Sedov solution in the non-relativistic 
limit. So in this paper, we will construct a dynamical model that is really capable of describing 
generic fireballs, no matter whether they are radiative or adiabatic, and no matter whether they 
are ultra-relativistic or non-relativistic. 



2 CONVENTIONAL DYNAMIC MODEL 

A differential equation has been proposed to depict the expansion of GRB remnants (Chiang & 
Dernier 1998; Piran 1998), 

*L = -Oj± (1) 

dm M ' y ' 

where to is the rest mass of the swept-up medium, 7 is the bulk Lorentz factor and M is the total 
mass in the co-moving frame, including internal energy U. Since thermal energy produced during 
the collisions is dE = (7— i)dmcP, usually we assume: dM = (l—e)dE/c 2 +dm = [(1— e)~/+e]dm, 
where e is defined as the fraction of the shock generated thermal energy (in the co-moving frame) 
that is radiated (Piran 1998). It is putative that equation (1) is correct in both ultra-relativistic 
and non-relativistic phase, for both radiative and adiabatic fireballs. However, after careful 
inspection, we find that during the non-relativistic phase of an adiabatic expansion, equation 
(1) could not give out a solution consistent with the Sedov results (Sedov 1969). 

2.1 Radiative case 

In the highly radiative case, e = 1, dM = dm, equation (1) reduces to, 

^L = ^~ 1 (2) 

dm M ej + to ' w 

where M e j is the mass ejected from the GRB central engine. Then an analytic solution is 
available (Blandford & McKee 1976; Piran 1998): 

(7-l)(7Q + l) _ ( m + M c j \ 2 

(7 + IX70-I) \m + M ej ) ' { } 

where 70 and too are initial values of 7 and m respectively. Usually we assume 70 ~ tj/2, 
itiq ~ M c j/r], where r\ = I*7o/(M e jC 2 ) and Eq is the total energy in the initial fireball (Waxman 
1997; Piran 1998). 

During the ultra-relativistic phase, )>1, M c j 3> m, equation (3) gives (7 + 1)to ~ M e j, 
or equivalently the familiar power-law 7 oc i?~ 3 , where R is the radius of the blast wave. In the 



later non-relativistic phase, 7~l,m> A^ej) we have: m 2 (3 2 = 4M 2 -, or (3 oc -R -3 , where (3 = v/c 
and v is the bulk velocity of the material. This is consistent with the late isothermal phase of 
the expansion of supernova remnants (SNRs) (Spitzer 1968). From these approximations, we 
believe that equation (2) is really correct for highly radiative fireballs. 

2.2 Adiabatic case 

In the adiabatic case, e = 0, dM = 7<im, equation (1) also has an analytic solution (Chiang & 
Dernier 1998) : 

M = [Ml + 2 7o M cj m + m 2 ] 1/2 , (4) 

m + 7 M ej 
7= M ■ (5) 

During the ultra-relativistic phase, 70-Wcj ^ fn S> Mcj/70, 7> 1, this solution can produce 
the familiar power-law 7 oc R~ s ' 2 , which is often quoted for an adiabatic blastwave decelerating 
in a uniform medium. In the non-relativistic limit (7 ~ 1, m S> 7oM c j), Chiang & Dernier 
(1998) have derived 7 ~ 1 + 7oM c j/m, so that they believe it also agrees with the Sedov solution 
(Lozinskaya 1992). However we find that their approximation is not accurate enough, because 
they have omitted some first-order infinitesimals of joM e j/m. The correct approximation could 
be obtained only by retaining all the first and second order infinitesimals, which in fact gives: 
7 ss 1 + (7oM e j/m) 2 /2, then we have (3 oc i?~ 3 . This is not consistent with the Sedov solution! 
We have also evaluated equation (1) numerically, the result is consistent with equations (4) and 
(5), all pointing to (3 oc -R~ 3 , not the relation of f3 oc R~ 3 ' 2 as declared popularly in the literature 
(Chiang & Dermer 1998; Piran 1998). 

This question is serious. First, it means that equation (1) is not a dependable model for 
non-radiative fireballs, although it can reproduce the major features in the ultra-relativistic 
phase. Second, the expansion of a realistic fireball is widely believed to be highly radiative 
at first, but after only several days, the expansion will become non-radiative (Sari, Piran & 
Narayan 1998; Dai, Huang & Lu 1998). So in the non-relativistic phase, the fireball is more 
likely to be adiabatic rather than highly radiative. However, it is just in this condition that the 



conventional model fails. So any calculation made according to equation (1) will lead to serious 
deviations in the light curves in the non-relativistic phase. 

3 OUR GENERIC MODEL 

Equation (1) is not consistent with the Sedov solution, we need to revise it. In the fixed frame, 
since the total kinetic energy of the fireball is Ek = (7 — l)(M c j + m)c 2 + (1 — e)^U (Panaitescu, 
Meszaros &; Rees 1998), and the radiated thermal energy is £7(7 — l)dmc 2 (Blandford & McKee 
1976), we have: 

d[(7 - l)(M cj + m)c 2 + (1 - e)-yU] = -€7(7 - l)dm c 2 . (6) 

For the item U, it is usually assumed: dU = (7 — \)dm c 2 (Panaitescu, Meszaros & Rees 1998). 
Equation (1) has been derived just in this way. However, the jump conditions (Blandford & 
McKee 1976) at the forward shock imply that U = (7 — l)mc 2 , so we suggest that the correct 
expression for dU should be: dU = d[(7 — l)mc 2 ] = (7 — l)dm c 2 + mc 2 d / y. Here we simply use 
U = (7 — l)mc 2 and substitute it into equation (6), then it is easy to get: 

*L = 1-zl . (7) 

dm M c j + em + 2(1 — 6)7771, 

We expect this equation should describe a generic fireball correctly. 

Indeed, in the highly radiative case (e = 1), equation (7) reduces to equation (2) exactly. 
While in the adiabatic case (e = 0), equation (7) reduces to : 

dj 7 2 — 1 



dm M c j + 27777, 

This equation has an analytic solution: 



(8) 



(7 - l)M cj c 2 + ( 7 2 - l)mc 2 = £ K o, (9) 

where Sko is the initial value of E^. In the ultra-relativistic phase (70-Mcj 3> m ^> Mej/7), we 
get the familiar relation of 7 oc i? -3 ' 2 . And in the non-relativistic phase (m 3> M e j), we get 
(3 oc i? -3 ' 2 as required by the Sedov solution. 



For any other e value between and 1, equation (7) describes the evolution of a partially 
radiative fireball. Unfortunately, we now could not find an exact analytic solution for equation 
(7). But in the non-relativistic phase, by assuming m 3> M e j, we still can get 771(7 — lp 2 ~ e '' 2 = 
const, that is: 

/?cxir 3 /( 2 - e ). (10) 

4 NUMERICAL RESULTS 

We have evaluated equation (7) numerically, bearing in mind that (Huang et al. 1998): 

dm = AttR nrripdR, (11) 



dR = /3c 7 ( 7 + ^J-f - l)dt, (12) 

where n is the number density of interstellar medium, m p is the mass of proton, and t is the 
time measured by an observer. We take Eq = 10 52 ergs, n = 1 cm -3 , M e j = 2 x 10 -5 Mq. 
Figures (1) - (4) illustrate the evolution of 7, v, R, and E-& respectively. In these figures, we 
have set e = (full lines), 0.5 (dotted lines), and 1 (dashed lines). It is clearly shown that our 
generic model overcomes the shortcoming of equation (1). 

For example, for highly radiative expansion, the dashed lines in these figures approximately 

satisfy 7 oc £~ 3 / 7 , R ex i 1 / 7 , 7 ex R~ 3 , Ek ex £~ 3 / 7 when 7 > 1, and v oc £~ 3 / 4 , R oc i 1 / 4 , 
v cc i?~ 3 , Ex x t~ 3 ' 4 when 7 ~ 1. While for adiabatic expansion, the full lines satisfy 7 x £~ 3 ' 8 , 
R x i 1 / 4 , 7 x i?~ 3 / 2 when 7 » 1, and satisfy v x £~ 3 / 5 , i? x £ 2 / 5 , v oc Z?~ 3 / 2 when 7 ~ 1. 

5 CONCLUSION AND DISCUSSION 

The conventional dynamic model is successful at describing highly radiative GRB remnants, 
however it has difficulty in reproducing the Sedov solution for adiabatic fireballs. This is com- 
pletely unnoticed in the literature. We have constructed a new generic model to overcome this 
shortcoming. Numerical evaluation has proved that our model is highly credible. We hope 



this work would remind researchers of the importance of the transition from ultra-relativistic to 
non-relativistic phase, which might occur as early as 10 6 - 10 s since the initial burst (Huang, 
Dai & Lu 1998). 

In the above analysis, for simplicity, we have assumed that e is a constant. But in realistic 
fireballs, e is expected to evolve from 1 to due to the changes in the relative importance of 
synchrotron-radiation- induced and expansion-induced loss of energy (Dai, Huang & Lu 1998). 
Assuming electrons in the co-moving frame carry a fraction £ e = 1 of the total thermal energy 
and that the magnetic energy density is a fraction £g = 0.01 of it, we re-evaluate equation (6) 
numerically. The results are plotted in Figs (1) - (4) with dash-dotted lines. We see from Fig. 
(4) that the evolution of e changes £k (t) dramatically. 

It is worth mentioning that SNRs evolve from non-radiative to radiative stage, but GRB 
remnants are just on the contrary. This is not surprising because GRB remnants radiate mainly 
through synchrotron radiation while SNRs lose energy due to excited ions. It is reasonable to 
deduce that at very late stages, when the cooling due to ions becomes important, GRB remnants 
may become highly radiative again, just in the same way that SNRs do. The transition may occur 
when the temperature drops to below ~ 10 6 K and the velocity is just several tens kilometer 
per second. This needs to be addressed in more detail. 

Another interesting problem is the possibility that HI supershells might be highly evolved 
GRB remnants (Loeb & Perna 1998; Efremov, Elmegreen & Hodge 1998). Our Figs (3) and (4) 
have shown that typical adiabatic GRB fireballs can evolve to R ~ 1 kpc at t ~ 10 6 - 10 7 yr, 
with v ~ 10 km/s, but highly radiative fireballs are obviously not powerful enough. To discuss 
this possibility in detail, we should pay attention to the possible adiabatic-to-radiative transition 
mentioned just above. 

We are very grateful to R. Wijers for his valuable comments and suggestions. This work 
was partly supported by the National Natural Science Foundation of China, grants 19773007 
and 19825109, and the National Climbing Project on Fundamental Researches. 
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Figure Captions 



Figure 1. Evolution of the bulk Lorentz factor 7. We take Eq = 10 52 ergs, n = 1 cm -3 , 
M e j = 2 x 10 -5 Mq. The full, dotted, and dashed lines correspond to e = (adiabatic), 0.5 
(partially radiative), and 1 (highly radiative) respectively. The dash-dotted line is plotted by 
allowing e to evolve with time (see Sect. (5) in the text). 



Figure 2. Evolution of the bulk velocity v. Parameters and line styles are the same as in Fig. 1. 

Figure 3. Evolution of the shock radius R. Parameters and line styles are the same as in Fig. 1. 

Figure 4. Evolution of the total kinetic energy £"k- Parameters and line styles are the same as 
in Fig. 1. 
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